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q", Pion production from a critical QCD phase 

o : 

O . N.G. Antoniou, Y.F. Contoyiannis, F.K. Diakonos, 
A.I. Karanikas and C.N. Ktorides a 

A theoretical scheme which relates multiparticle states generated in ultrarelativistic 
nuclear collisions to a QCD phase transition is considered in the framework of the uni- 
versality class provided by the 3-D Ising model. Two different evolution scenarios for the 

■^J- ! QGP system are examined. The statistical mechanics of the critical state is accounted 
for in terms of (critical) cluster formation consistent with suitably cast effective action 

(SI ' functional, one for each considered type of expansion. Fractal properties associated with 
these clusters, characterizing the density fluctuations near the QCD critical point, are 
determined. Monte-Carlo simulations are employed to generate 'events', pertaining to 
the total system, which correspond to signals associated with unconventional sources of 
Ph| pion production. 

1. Introduction 

It is widely speculated that ultra-relativistic nucleus-nucleus (A+A) collisions offer 
^ ■ an experimental methodology by which a transition from ordinary (confined) matter to a 
quark-gluon plasma (QGP) state can be transiently attained. In view of such a prospective 
it becomes imperative to identify specific experimental signals in the outcome of these 
collisions whose origin can be unmistakenly ascribed to a critical behaviour. 

From a phenomenological point of view we can list three basic topics as the most 
relevant to such an investigation: 

(a) How to pick up patterns in the multihadron (pion) states generated in A+A collisions 
from the central rapidity region and extract pertinent information from them. 

(b) What is one's conception of the evolution picture taking place between the (pre- 
sumed) formation of the quark-gluon plasma system and hadronization. 

(c) In what way does the underlying microscopic theory, i.e., QCD, leave its imprint on 
the production patterns of hadrons. 

There is a variety of attitudes that can be taken with respect to the above issues. Even 
though the present work focuses upon the last one, it is important to specify the adopted 
viewpoints on the other two as well. 

In regard to (a), our premise draws its inspiration from the Bialas-Peschanski |I| ap- 
proach to data analysis for large multiplicity events. The emphasis, here, lies on the 
possible presence of non-statistical deviations which can be uncovered from experimen- 
tally determined probability distributions. The idea is to look for big bursts from small 
regions (cells) of the phase space which signal unusual fluctuations of density, an occurence 
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that has been termed as search for intermittency. Such intermittent spiky events induce 
a dependence of the scaled moments of experimentally defined probability distribution 
P(pi, ■■■,p n ) on the size SQ of each cell, according to some power-law. Searches for inter- 
mittency in multiparticle states resulting from A+A collisions have so far revealed only 
fluctuations due to Bose-Einstein interference. However, the question that still remains 
open is whether intermittency could be directly ascribed to a quark-hadron phase transi- 
tion of second order. This is precisely the attitude that underlines our present approach, 
i.e. our intention is to propose a scheme according to which intermittency patterns in 
a multiparticle hadronic state are intimately linked with fluctuations which make their 
presence near the QCD critical point. 

Concerning, next, the space-time evolution of the (plasma) system produced in the cen- 
tral rapidity region we shall adopt, for the most part, the picture promoted by Bjorken 
[0,0. The basic idea is that, following the collision, a "central plateau" is formed com- 
prised of a thin slab which contains quarks and gluons in (local) thermal equilibrium. The 
expansion is envisioned to be smooth enough so that conditions of thermal equilibrium 
persist on a locus which is monitored by a collection of inertial observers with relative 
velocities. A proper time (r) dependence of all relevant thermodynamic quantities that 
describe the "fluid" emerges so that one writes e (r), p (r), j3 (r), respectively for energy 
and pressure densities and for inverse temperature. Dynamical evolution is thereby de- 
scribed by a sequence of a hyperboloidal, spacelike surfaces whose profile is specified by 
a hyperbolic section in the t — z plane (z is the direction along which the collision takes 
place) and flat geometry in the transverse directions. This situation persists for a time 
scale of ~ 10-20 fm/c (~ /? c ^4 3^ at which point formation of primordial hadrons (mostly 
pions) occurs. As an alternative possibility to the above we shall consider, in Section 5, a 
spherically symmetric evolution of a centrally formed fireball of quarks and gluons. Our 
numerical studies will cover both scenarios. 

The final issue, towards which the present paper basically addresses itself, concerns 
the interpolation between the underlying microscopic theory and observed patterns in the 
multihadron states generated by A+A collisions. In particular, we shall focus our efforts 
on determining possible imprints left by long range fluctuations, present in a critical 
quark-gluon system, on the produced multiparticle states. Evidence is, in fact, mounting 
from ongoing theoretical activity based on microscopic considerations which simulate QCD 
through specific models (e.g. that of Nambu and Jona-Lasinio), which lends overwhelming 
support to the existence of a critical point in the thermodynamical plane of temperature 
vs. chemical potential @j7j. The latter is the remnant of a tricritical point occuring in 
the phase diagram of the ideal case m u = = and from which a second order chiral 
phase transition curve originates belonging to the universality class of the 0(4) a- model 
[||] . For the realistic case where the up-down current quark masses are different from zero 
it marks the termination of a first order phase transition curve 0. Such a point signals 
the onset of long range correlations and the relevant issue is to identify its universality 
class. It is widely believed that the latter is represented by the 3-D Ising model 
Under this assumption we shall adopt a strategy according to which the corresponding 
effective action T [a] is: a) accomodated on hyperboloidal surfaces of the type specified by 
Bjorken's picture for the space-time expansion of the plasma or b) adjusted to a spherical 
expansion scenario. Given that Y [a] is a functional of a classical field, having resulted 
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from an integration over microscopic degrees of freedom, we shall propose that it plays 
the role of the free energy associated with sigma condensates. 

From such a premise, the bulk of our effort will be to establish direct connections 
between the critical QGP system and observed patterns in hadron production. It is 
our aim to make specific predictions relating to the extensive experimental programme on 
heavy ions now in progress (SPS, RHIC, LHC), some of which have already been reported 
in previous, brief, expositions ||. Our results will be summarized in an assortment 
of Monte-Carlo simulations stemming from the advocated theoretical approach. Given 
that intermittency effects in hadron production patterns are foremost in our minds, the 
significance of the presented output should be directly relevant to the goal of identifying 
clear experimental signals for the existence of a critical point in quark matter. 

As can be surmised from the above remarks, our proposed scheme relies centrally on 
the idea that the observed patterns in hadron production have a component whose source 
is the large scale dynamics induced by the underlying microscopic theory at the critical 
point. At the same time, it constitutes an effort which stays close to an equilibrium 
frame of description^, in juxtaposition with alternative approaches which place primary 
emphasis on the non-equlibrium aspects of the problem, see, e.g. | |10[ |. 

Clearly, we cannot altogether avoid the out-of-equilibrium dynamical component enter- 
ing the transition from the (critical) QGP to the hadronic state of the considered physical 
problem. In our scheme this matter enters in connection with the critical slowing down 
expected to occur between the critical, T c , and freeze-out, Tf, temperatures - the latter 
marking the onset of the emergence of hadrons. Adopting a working assumption of the 
form T c ^ Tf, owing to the fact that the specific heat diverges at the critical point, allows 
us to remain within the equilibrium frame of description. On the other hand, the (proper) 
time interval between the aformentioned two temperatures is not necessarily negligible. 
Clearly, dynamical matters related to critical slowing down effects present an interest of 
their own. We shall not, in this paper enter the discussion of issues related to dynamical 
scaling behavior. 

Our paper is organized as follows. Section 2 is devoted to the presentation of the field 
theoretical framework within which we shall pursue the study of the critical system, as- 
suming Bjorken's evolution picture for the QGP. In Section 3 we develop a Statistical 
Mechanical frame of description for the critical field system which relies on cluster for- 
mation. We proceed, in Section 4, to discuss the statistical content of the global system, 
setting the stage for the development of a Monte-Carlo generator for critical events that 
will be taken up in Section 6. In the interim, Section 5, we shall display the basic battery 
of formulas pertaining to the spherical evolution scenario for the QGP system. Section 
6 addresses the quintessential issue regarding the generation of events with critical fluc- 
tuations as well as their identification in the profile of the produced pions which can be 
traced to a critical point of second order. In section 7 we present our numerical results 
by performing an event-by-event analysis in a large set of critical events and showing how 
our dedection algorithm of the critical fluctuations works in practice. Finally, section 8 
contains concluding remarks. 

1 We shall not, in this paper, be preoccupied with precise dynamical accounts underlying the time devel- 
opment of the QGP system from its initial formation to its critical stage. Our basic concern is, simply, 
to justify our choices for the geometrical profile of the QGP state as the critical point is approached. 
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2. Adjustment of the 3-D Ising-model to Bjorken's evolution scenario 



Our basic theoretical input is provided by the effective action of the 3-D Ising model 
r[er], which, following the analysis of [TT|, adequately represents QCD in the vicinity of 
the second order critical point T = T c , fi = fi c (as the endpoint of a first-order phase 
transition line) resulting by taking into account the finite mass for up and down quarks. 
Thus, provided driving parameters of the model (e.g., temperature, coupling constant, 
etc.) are on or near their critical values, r c [a] can be looked upon as the free energy for 
the critical system. 

According to our viewpoint the aformentioned free energy incorporates the sigma field 
condensates, a ~ (qq), in thermal equilibrium near the critical temperature T w T c . As 
already mentioned in the introductory section, for the most part of this work we shall feel 
phenomenologically justified to keep our considerations on a strictly equilibrium level by 
taking into account that the freeze-out temperature Tf, at which one meets the asymptotic 
states of the produced hadronic system, is very close to the critical temperature (T c ~ Tf). 
This means that we shall view the second order phase as a direct interpolation between 
the condensates and the asymptotic hadronic states. 

Under conditions of strict local equilibrium the thermodynamics of the condensates is 
described, in the ordered phase (T < T c ), by the free energy which we express, at T = T c , 
as follows 



r>] = t- 1 J <?x[\{VoY + ^(r-v; 



5+n 



(1) 



where 5 is the isothermal critical exponent and g is a dimensionless coupling. The field 
a has a dimension a ~ (length) _1 . The parameters (g,S) are universal and specify the 



For the Ising 3-D universality class we have S 



and g « 1.5 — 2.5 as obtained in 



T c [a] in a manner that is explicitly 



critical equation of state: ~ go 
due to the smallness of the anomalous dimension 

[00- 

Our next concern is to cast the effective action 
adapted to Bjorken's inside-outside cascade picture for the plasma expansion. For this 
purpose we introduce rapidity and proper time coordinates (£, r) so that the longitudinal 
space element, corresponding to a local observer in a comoving system, becomes dx\\ = 
rd^. With this choice and for the purpose of describing the clusters formed at T = T c , the 
"longitudinal" integration occurs along the critical hyperbola r = r c . We thereby obtain 
the following expression for the effective action which furnishes an appropriate account 
for the system as the critical temperature is reached: 



r c [c] = 77- / d£ cosh £ / d 2 x ± 



2 cosh £ 




(V ± a) 2 + <7^Cl(T c -V; 



5+1 



where A is the rapidity size, 5*j_ is the transverse area of the system and Ca = j-- 

One notices that, through (Q), the system has been accomodated on a space-like hyper- 
boloidal surface which is flat in the transverse directions x, y and intersects the t — z plane 
as a hyperbola (z is the collision direction). The monitoring of the field system in our case 
calls for a collection of local inertial observers who must put their data together in order 
to construct its full description. This occurence is ideally suited to the Bialas-Peschanski 
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scenario concerning the detection of intermittency patterns in the multiparticle hadronic 
production. Within this picture it turns out that the total rapidity size A in eq. (fj) should 
be replaced by a narrow strip A£ in rapidity around a local observer. Furthermore, we 
are free to choose A£ < 1 and make use of the approximation cosh£ « 1 which is valid 
within 10%. 

As already stated our central goal is, assuming that T c « Tf, to formulate correspon- 
dences between field quantities associated with the critical system, on the one hand and 
hadronic observables on the other. On the field theory side the relevant issue concerns 
the statistical account of the macroscopic system comprised of condensates. One needs 
to define a statistical weight factor appropriate for the equilibrium situation in hand. A 
convenient way in terms of which one might proceed to accomplish this task is to employ 
a coherent state analysis (see, e.g., |L4||) which, by definition, pertains to collective fields. 
Let us recall that coherent states are eigenstates of field operators, being a superposi- 
tion of Fock states of arbitrary particle content. One writes (the tilde on top signifies an 
operator) 

a(x)\a > = a(x)\a >, (3) 

where the \a > denote coherent (classical, in effect) states of the field a. Our basic 
assumption is that the density matrix associated with the partition function is diagonal 
in the coherent-state representation. 

In turn, the average multiplicity < n > of sigmas is given by 

<n>= J V[a] J v d 2 x ± d£a 2 (x ± ,0 e"^ 1 / J V[a] e - r ^ ] (4) 

The above formulas will be centrally relied upon in our subsequent analysis. They, in 
fact, serve as a bridge between an original microscopic description - which has been given 
a macroscopic, collective content by integrating out a huge number of degrees of freedom 
- and the statistical treatmeant of fluctuations in the resulting critical, collective system 

m. 



3. Statistical Mechanics of the critical system: cluster formation 

To the extent that we shall be restricting ourselves to a strictly equilibrium mode of 
description, the strategy by which we propose to study the critical system is through 
cluster formation. The physical picture we attach to (critical) clusters is that they rep- 
resent distributions of primordial hadrons (massless sigmas) whose growth is driven by 
the transition from the QGP to the hadronic state of matter. From a formal, mathemati- 
cal, viewpoint such clusters serve as entities which interpolate between critical exponents 
characterizing the plasma system on the one hand and fractal dimensions associated with 
the hadronic (pion) emission patterns on the other. 

The 3-D configurations representing the clusters will be assembled as follows. First, 
their one-dimensional (longitudinal) profile along the rapidity axis will be considered, 
given that this is where the bulk of the dynamical effects is encountered. Next, we shall 
proceed to investigate their transverse fractal properties by freezing the rapidity variable. 
The cartesian synthesis of a given cluster, provides a well defined picture of a cluster's 
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constitution. The latter, certainly adheres to the cylindrical (with respect to the rapidity 
axis) and Lorentz (save at the edges) symmetries which characterize the overall system. 

The projection onto the rapidity direction is furnished by the 1-D effective action func- 
tional: 



a 




a 



S+1 
2 



(5) 



where A£ is the rapidity strip around the local observer at £ = and R± the transverse 
radius of the system at T = T c , while a denotes a dimensionless cr-field (a = T^a). 



For notational ease we make the identifications 
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(i) 



PcT c 



a (1) 



which recast Eq.([|) into the following generic form 

2 



cr 



si" 



AS 



HI) ^ 



a 



(6) 



(7) 



The local system is considered as open, i.e. it communicates with the total system, so 
that no boundary conditions are imposed on it. In turn, this implies that we are free 
to consider analytic continuation of (|7]), with respect to A£, which will furnish scaling 
properties of extensive thermodynamical quantities registered by the local (sub) system 
as its size is being varied. We emphasize that our specification £ = has been made 
strictly for notational ease; any local oberver, located at £ = £j is equally acceptable, due 
to Lorentz invariance, provided that the said observer is not too close to the boundaries 
of the total system. 



For heavy nuclei, 



(3cTc 



^> 1, the evaluation of the partition function of the subsystem 



V[a]e 



(8) 



is dominated by saddle point configurations. In what follows we drop, for simplicity of 
notation, the hat used to indicate dimensionless quantities. 

There emerge instanton-like solutions of the classical equations of motion (ID 

a -{6+ 1)^V = 
which are classified according to total energy 



(9) 



gi\\ s+1 



(10) 



and size £o- 

Now, for a given solution it follows, in a straight forward manner, that a suppression 
term e~ EA ^ factorizes in the expression for the partition function, hence the dominant con- 
tributions correspond to configurations with vanishing energy. Under these circumstances 
one obtains the following analytic expression for the instanton-like solution 
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As a function of the varying size A£ the effective action takes the form 
2 9 tiR 2 

(3 C T C ~JA£- 



r c (A£ ; £ 



-C\\ H£-£o|] m ^ 



(12) 



In order to suppress large contributions to the free energy, hence practically zero con- 
tributions to the partition function, we impose the restriction that the "instanton" size 
be much larger than the size of the local system, i.e. A£ <C £o- As a result, we obtain 
approximately constant solutions in the domain A£, of the form 



a 



V2 



2 
4-1 



(13) 



The computation of the partition function is now reduced to an integration over the 
"instanton" size. We let physical insight guide the relevant calculation by thinking as 
follows. Consider an extensive variable M associated with the field configuration of the 
critical system in the region monitored by our local observer, including its analytic contin- 
uation. Specifically, we choose the quantity M = J^[a(x)] 2 dx which furnishes, according 
to eq.(f|), the multiplicity < n(A£) >=< J^[a(x)] 2 dx > of sigmas within the domain of 
extent A£. Moreover, we introduce the concept of a cluster of radius A£ as an object with 
geometrical properties built up through the statistical average over configurations corre- 
sponding to values of M greater than a minimum value \x (for the case at hand \i — 1). 
Referring to (f|) and flT2"| ) this average is determined as follows 



< 



[a(x)] 2 dx > 



where d = 2g{ 1] g^ 1 A s+1 and A = [g ( 2 L, /2(5 - 1)T^. 

It can now be explicitly shown |T7[ that, for G\ 3> 1, three characteristic regions can 
be clearly identified with respect to the behaviour of the integral on the right hand side 
of (0): " 




4 + 1 
4-1 



a- 1 - (Co - ao 



xexp 



AO" 



4+3 



4+3 1 
" 4-1 



(14) 



A£< A d 
A d < A£ < A u 
A£ > A u 



< n(A£) >~ const. 
<n(AO >~ (AO^ 

< n(Af) >~ (AO^ 



(15) 



where 



4+1 



A d = A-^Gf- 1 ^ 



A u = G{ 



4-1 
4+3 



(16) 



serve as the lower (A d ) and the upper radius (A u ) of a cluster centered around £ = 0. 
There follows a fractal structure for the critical clusters with mass dimension 

5-1 



d 



(i) 



5 + 1 



(17) 
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Turning our attention to the two-dimensional, transverse profile of the critical system 
we encounter the effective action functional 



rf[<7] = C A A J d 2 x ± \l(V ± a ± y + g (a 2 ) 



6+1 

2 



;is) 
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where _L, as subscript to the integral, serves to denote transverse size and A furnishes the 
total rapidity. 

The form of the above expression is a two-dimensional analogue of (H) with correspond- 
ing parameters 

g{ 2) = C A A ; gi 2) = g (19) 
The counterpart of eq.flToT) now reads 

R <^ Rd ; < n±(R) >~ const. 

R d < R < R u ; < n±(R) >~ R 2 ^ 1 

R->R U ; < n ± (R) >~ R 2 ~^ (20) 

where 

R d = GfA^ ; R U = (3 C G^ (21) 

with G 2 = 2irg ( ? ) gi 2) A s + 1 (*±2) and A 2 = ((#/4)(<5 - 1) 2 (5 + In analogy 

with eqs.flloD we conclude from eqs.(^) the formation of fractal clusters in transverse 
space with mass dimension dp = ^pp- in the range Rd <C R <C R u - The 3-D cluster is 
constructed following our original considerations according to which this object is roughly 
confined within a cylinder with radius R u in transverse space and size 2A U in rapidity. 
This realization has the cartesian product form (rapidity (g) transverse space) and it is 
consistent with longitudinal expansion of the original system. 

Referring to eq.([|), with T c [a] given by fl5|), we estimate the average multiplicity for 
the 3-dimensional cluster to be |T3[ : 

r ( 3 ) y f^l 

< » >-= (vi) (22) 

with 7 = 27ri22A u and V a = P*y/2gU2- 

4. Multi-cluster description of the global system 

Assuming that the entire critical system is a cylinder with rapidity size A and transverse 
radius R± we can calculate the mean number of 'cylindrical' clusters with rapidity size 
2A U and transverse radius R u contained in the global system. Denoting this number by 
N we determine N — N\\ x N±, where N\\ denotes the number of clusters in rapidity and 
N± in transverse space, respectively. Using the fact that k = ^±1 m 3 for the Ising 3-D 



universality class [|12l , |i3|| we estimate, from the previous analysis, these numbers as: 
11 2A U ^R ± \2J ' 
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It follows that the 3-D critical system appears as an almost cylindrically symmetric^ 
arrangement of, in the mean, N clusters. Each cluster (of size 2A U in rapidity and R u 
transverse radius) has a fractal mass dimension dp decomposed as the cartesian product 
dp = dp + dp of two fractals in the corresponding subspaces. We also observe that the 
size of the clusters in rapidity is relatively small, A u = ^ ("T*)' justifying a posteriori 
the approximation cosh£ « 1 made in eq.(Q). 

The complete analysis calls for the power law underlying the behavior of the density- 
density correlation function p(£, 0) =< p(£)p(0) > in the scaling region. Differentiating 
the mean multiplicity, given by eq.(22), with respect to rapidity we are led |j| to 



p(£,o) = ( 



27ngC A 



) 1/3 &" 



,£(1/2) 

r(i/6) 



-1/3 1 X -L 

1 /5c 



-2/3 



(24) 



As already mentioned, the fractal dimension emerging from the power-laws in rapidity 
space is given by d!p) = ^ and is valid in the scaling region 8 (= A d ) <C |£| < A u . 

In the presence of Nn non-overlapping sources, located along the rapidity direction at 
random points £i, ...,^ N , the corresponding partition function reads: 



iVn-l 



e n 



r/t- 1 



Pi 



I d^<a 2 >] 



(25) 



which leads to the following distribution for the given multiplicity Nf 



' A + 26 



i 



n 



(26) 



< f . One 



where 7 = |T , with the assumed ordering being — A < ^ < ... < ^ 
readily notes that the largeness of 7 allows us to regard the chain of clusters along the 
rapidity direction as a collection of practically non-interacting objects. This occurence 
will greatly facilitate, later on, our numerical simulations. 

Unlike the rapidity case, fluctuations associated with transverse space are not directly 
observable in heavy ion experiments. What actually is observed is the pion distribution in 
the transverse momentum space. The need, therefore, arises to transform the geometrical 
picture valid for the transverse space to the corresponding picture in transverse momentum 
space. The two spaces are connected at the level of the density-density correlation through 
a Fourier transform: 



d 2 x±e 



p(x ± ,0) 



(27) 



where p(x±,0) is the density-density correlation for the sigmas in the transverse space 
projection of a critical cluster with center at x' ± = and p(pp) the corresponding density- 
density correlation in the transverse momentum space. The latter is, by definition, related 
to the gradient of the mean multiplicity < n±(R) >, see eqs.(|2~0D, with respect to R. 



2 We say 'almost' because the individual clusters are not supposed to behave as rigid structures; they can 
deform, due to the presence of neighbouring clusters, both in transverse and in rapidity directions. 
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Setting R = \x±\ we obtain 



Performing the Fourier integration, according to (|27|) , we arrive at the following power-law 



for the density-density correlation in the transverse momentum space: 

2(k-1) 

PKPt) ~\Pt\ k (29) 

In analogy to the space description of the 3-D system, we now consider a collection of 
cylindrical clusters in the transverse momentum (g) rapidity space. The number of clusters 
in the transverse momentum projection is Nj_ given by eqs.(|2"3"|). Finally, the multiplicity 
of hadrons within one cluster is given by eq. (|22|) . 

In order to complete the description of the critical system in the transverse momentum 
we need an additional phenomenological input. We adopt the following assumptions: 

• The centers of the clusters in the transverse momentum space px,% are random 
variables distributed according to an exponential law: 

2 2 PT 

f(PT) = — ~i e <PT> 30 

7T < Pt > 

• The net transverse momentum flow through the clusters is suppressed: 

X>.< = ° ( 31 ) 

1=1 

where < px > is the mean transverse momentum which can be approximated to: < pr >~ 
2T C = 280 MeV. The size |Ap_j_| of each cluster in the transverse momentum space is 
determined by the minimum length scale Rd in the transverse configuration space and is 



given by the relation [15 



|ApJ = r 3/8 (^) V4 (32) 

Having exhibited the portrait of the critical system according to the longitudinal evo- 
lution scheme, we can now turn our attention to the spherical mode of expansion. This 
we shall do in the following section through a display of a minimum of formalism. 

5. Spherical evolution 

In the present section we shall consider a different evolution scenario according to which 
the geometry of the A+A collision is spherical. More specifically, we consider the situation 
where, following the collision, a fireball of spherical size is formed which expands radially 
outwards. As a result, the sigmas are now distributed inside a spherical cluster. 
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The starting point of our analysis is, once again, the 3-D effective action (1). Following 
a similar procedure with the one applied for the cylindrical geometry we find that the 
average multiplicity < N a > is given by 

< N a >= Z- 1 JvWe-^ip- 1 J d 3 ra 2 (r)) (33) 

where the insertion of the factor (3~ l has been made on dimensional grounds since dim[a 2 ] = 
(length)~ 2 and < N a > must be dimensionless. 

The partition function at T = T c , with the scalar particles (a) in a spherical volume 
V{< (4/3)7i\R^), now has the form 

Z(N a , V, T c ) = N- 1 ' 2 expH-^r 1 ^) (34) 
where V = {2g)^ fil 



The upper limit R u of the cluster can be computed directly from [17]. For the D = 3 



case the anomalous dimension r\ of the sigma-field is involved in the resulting expression. 



i 



We obtain R u = G 3 q with q = 3- ^±§^ and G 3 = f^^gf ] Q- F ° r small r] (77 < 1) 

we find that R u — > 00. Accordingly, we can consider the whole system as being comprised 
of only one cluster with radius R u . The average multiplicity < N a > in this cluster, as a 
function of the volume V, follows from fl34j) and is given by 

The above equation implies that there is a minimal volume scale, V = V Q , where self- 
similarity breaks down. 

On general grounds, the fractal dimension in i? 3 -space is given by dp = 3^- if we 
neglect 77. Taking rj into account (~ 0.034 [|13|]) we obtain dp = 1.98. Thus if we were 
to project this fractal onto rapidity and transverse space, respectively, we would not find 
corresponding fractal configurations since the dimensions of the spaces upon which we 
project are smaller or nearly equal to dp [ p~8j| . 

To complete our considerations in configuration space we give the expression for the 
density-density correlation function which obeys the following characteristic power-law 



*r 3 



< P(r)p(0) >= ^M(^)-«p—r-« (36) 



Turning our attention to momentum space, the pattern is again a spherical cluster with 
center at p — 0. The corresponding radius is determined by the lower limit Rd of the 
cluster size in configuration space, i.e. | Ap| ~ R^ . On the other hand, Rd is determined 
by the minimal volume V Q = (4/3)irR d . It follows that: 
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and for g pa 2, we estimate |Ap| pa 150 MeV. The fractal dimension in momentum space 
is dp = 3 — dp = 3(1 — s= ^) and, the momentum distribution of sigmas inside the cluster 
is determined through the density- density correlation: 

< P(b1)p(0) >~ bT 3 ^ (38) 

which results from the Fourier transform of eq . (|36|) . 

This concludes our statistical mechanical account for the critical fluctuations in terms 
of cluster formation, for each of the two alternative evolution schemes. We shall proceed, 
in the next section, to discuss corresponding numerical studies, pertaining to a critical 
component for pion production in A+A collisions through simulations of a Monte-Carlo 
type. Our expectation is that the 'theoretical generation' of events whose origin lies 
exclusively in a critical state reached by the QGP system, will aid the search for respective 
signals in produced multiparticle pion states. 

6. Monte-Carlo simulations 

Following the considerations in the previous sections it becomes obvious that the for- 
mation of critical clusters is expected to induce direct phenomenological consequences, 
in the multiplicity distribution of the produced sigmas, in the final states of a A + A 
collision. In the first part of the present section we shall describe the development of a 
Monte-Carlo programme for the purpose of simulating the production of "critical pions" 
in such collision processes, coming from the decay of light sigmas. The reader who is not 
interested for the technical details of our approach can skip the first subsection and go 
directly to the next one where we discuss a proposed algorithm for detecting the critical 
fluctuations in experimental data. 

6.1. CMC Generator for critical events 

In the context of the QCD phase diagram in the temperature/chemical potential plane, 
a two step procedure is required for the study of the critical pion spectrum formation. 
First, we need to produce the clusters corresponding to the cx-condensates. According to 
our theoretical discussion, the geometry of the cr-clusters is determined by the evolution 
pattern of the centrally formed fluid of quarks and gluons. In the case of longitudinal 
evolution we have a cylindrical arrangement of several clusters in the transverse momen- 
tum rapidity space (cf. sections 3 and 4) while in a spherically symmetric evolution a 
single cluster in momentum space is formed (cf. section 5). 

Second, we let the sigmas in the condensates decay into pions with a branching ratio 
1:2 for neutral to charged pion production. The crucial parameter for determining the 
spectrum of the produced pions is provided by the mass of the decaying sigmas. As the 
sigma mass is not a constant but, in fact, evolves during the freeze-out process, we treat 
it as a varying parameter in our numerical considerations. 

We recall that the number of clusters in the case of cylindrical symmetry is determined 
by the size of the critical region while the multiplicity of cr's within each cluster depends 
on the couplings occuring in the free energy of the system at the critical point. Once the 
sigma mass is assigned a specific value, we use two-body phase space kinematics in order 
to determine the transverse momenta and the rapidities of the pions. For the spherically 
symmetric case, on the other hand, a single critical cluster is formed around the CM 
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of the system (p = 0) while the multiplicity within the cluster depends, once again, on 
the couplings of the "critical" free energy. In what follows we shall conduct numerical 
investigations, using Monte-Carlo methodology, for both the cylindrical and the spherical 
configurations of the critical clusters. It could very well be that the geometry of the 
evolving critical system is neither perfectly cylindrical nor spherical but consistent with 
a picture interpolating between these two extreme regimes. 

Let us consider the implementation of the first step (simulation of the sigma-condensates) 
for a longitudinally expanding critical system (cylindrical geometry). The input parame- 
ters determining the sigma-distribution are: the size R± of the system in the transverse 
direction, the total rapidity interval A, the proper time scale r, the critical temperature 
T c , the coupling g in the 3-D Ising effective free energy and the isothermal critical ex- 
ponent 5. Once we fix the values of the input parameters we can calculate the size of a 
cylindrical cluster of sigmas, meaning its transverse momentum radius |Apj_|, given by 
eq. (p2[) , and its size in rapidity 2A U (see eq.flT6|)). The number of sigma clusters in the 
transverse (momentum or real) space projection of the critical system, as well as the cor- 
responding number of clusters in rapidity space are then given through eqs.(|23|). Finally, 
the multiplicity of sigmas within each cluster can be obtained using eq. (E£3) . 

Next, we construct the phase space projections of the sigmas in the critical event. The 
determination of the rapidity space profile proceeds as follows: We produce a configuration 
of N\\ centers for these clusters treating them as random variables, uniformly distributed in 
[0, A] according to eq.( p6|) . Then we generate around each center the sigmas, distributed 
with a density- density correlation following the power-law ([15]). 

This stage of the algorithm requires the construction of a random fractal with a given 
fractal dimension dp. In order to perform this construction we have to use test-functions 
of the form P(x) ~ x~ 1 ~ dp as a generalization of the Levy distributions for arbitrary 
embedding dimension |I9|. We produce randomly the rapidities corresponding to one 
cluster using the appropriate test-function. In general, the resulting set does not have the 
correct center and size. Due to self-similarity the cluster size is obtained through a simple 
scaling while, with a suitable translation, the cluster is placed around the corresponding 
center. Coalescence of neighbouring clusters is allowed, so the actual multiplicity within 
each cluster depends on the position of the centers of its neighbours. Therefore, the cluster 
arrangement in the rapidity space determines the total multiplicity N a of sigmas in the 
event. 

Turning our attention to the transverse momentum space we arrange the centers of 
the iVj_ clusters so that they are distributed according to the exponential law (|30l), under 
the constraint that the total transverse momentum of the centers vanishes. To each 
such cluster correspond ^r- sigma-particles distributed according to the power-law (pOj) . 
The construction algorithm for the transverse momentum clusters is analogous to the 
algorithm used for the construction in rapidity space. The test-function has to be adjusted 
to produce the random fractal with the correct dimension. 

We, therefore, end up with two sets of phase space variables. A set Sp with N a trans- 
verse momentum variables { (p T , Xl , Pt, Vi ) , (pi> a , Pt, V2 ),•••, {PT,x Na , VT )VNa ) } and a set 
with N a rapidity variables £ 2 > ■■•> 6v CT }- The sigma-content of a "critical" event, re- 
specting the cartesian product form of the underlying cylindrical geometry, is then realized 
as any possible one-to-one pairing of the elements of Sp with the elements of S%. 
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For a spherically symmetric expansion the size of the critical cluster is represented by 
a single parameter R corresponding to the radius of the (growing) system in the 3-D 
space. This makes the case in hand much easier to handle. First we specify the size 
| Ap\ of the (single) critical sigma cluster in momentum space, cf. eq.([J7]) and then using 
eq . (|35|) we determine the total sigma multiplicity N a within this cluster. The generation 
of the momenta of the sigma-particles is then based on the observation that for the 
3-D system the density-density correlation function in momentum space coincides with 
the step distribution of a Levy flight leading to a fractal set with the same fractal 
dimension (d F ~ 1). Thus the algorithm used for the generation of the a-momenta is 
given through the following equations 



Px,i = Px,i-i + |Ap]nyl - r|cos(27rr 3 ) 
p y ,i = p y ,i-i + \ Ap\n\Jl - r\ sin(27rr 3 ) 

Pz,i = Pz,i-i + \Ap\nr2 (39) 



where r 1 is distributed according to the density-density correlation function (|38|), r 3 is 
random variable uniformly distributed in [0, 1] while r2 is uniformly distributed in [—1,1]. 
A final scaling and translation, as in the case of the cylindrical clusters, is needed in order 
to embed the cluster in a region with radius \Ap\ around p = in momentum space. 

In order to establish contact with the observables in an experiment with colliding rela- 
tivistic heavy ions we must take into account the decay of sigmas into pions. As already 
discussed this constitutes a straight forward step for a given value of the mass of the de- 
caying sigmas. If this mass is far from the pion production threshold 277?.^ the kinematics 
of the pions is strongly influenced by the sigma-mass and the fluctuation pattern of the 
pions is disordered relative to the spectrum of the sigmas. On the other hand, the sigma 
mass is generated dynamically during the freeze-out phase and, if the decay rate is much 
larger than the expansion rate, there is a good possibility that the sigmas decay imme- 
diately after their effective mass has reached the pion production threshold fll5| . If this 



is the case then the fluctuation patterns describing the sigma momenta are transfered, 
almost unchanged, to the momenta of the produced pions. Here however, we will follow a 
more general mechanism for the decay of the sigmas allowing us also to describe a more 
conservative scenario when the sigmas decay well above the two pion threshold. Specifi- 
cally, in order to simulate the variation of the sigma mass during the freeze-out process, 
we introduce a probability distribution P(M a ) determining the number of sigmas with 
mass M a which decay into pions. We assume that P(M a ) is well described by a Gaussian 
of mean value M a and variance 5m a . 

The generation of the critical events is therefore completed in the last stage of our 
Monte-Carlo code (CMC) when, treating the sigma mass as a random variable with 
propability density P(M CT ), we produce the final charged pions through the decay of the 
critical sigmas. The quantum mechanical amplitude describing the probability of the 
sigma-decay into a 7r + -7r~ pair is supposed to be independent of the sigma-momentum. 

6.2. Reconstruction of the critical cr-sector 

Using the CMC algorithm we produce a large number of critical events and we study 
their phenomenology. Our main interest is to reveal the critical fluctuations from the 
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momenta of the final observed charged pions. For this purpose we perform an event- 
by-event intermittency analysis in the factorial moments of the negative pion momenta 
looking for a region with self-similar fluctuations. If the mass parameter M a is less than 
1 MeV above the two pion threshold the critical fluctuations can be directly observed 
in the momenta of the negative pions. Increasing M a the effect gets more and more 
suppressed leading, relatively soon (M a ~ 10 MeV + 2m,,-), to an almost dissapearance 
of the critical fluctuations (see Fig. 5, open up triangles). This behaviour remains even 
if the variance 5m a vanishes. If we increase Sm a for a given value of M a we find, in 
general, a suppression of the critical region with rate depending on the value of M a . If 
M a < 1 MeV + 2m n the rate of suppression is very small and the critical fluctuations, for 
realistic values of 5m a , practically persist. In order to get back the critical behaviour one 
has, therefore, to reconstruct the critical sigmas based on the momenta of the observed 
charged pions. We propose the following reconstruction algorithm: 

• For a given event, using as input the momenta of the charged pions, form the 
invariant mass m^. +7r _ = {p + + pJ) 2 for all possible pairs of 7r + -7r~ {p + (p_) is the 
4- momentum of 7r + (7r~)). 

• Determine the probability distribution of the quantity m 2 +n -. Look for a peak 
in the distribution. Form a new set of charged pions picking up the 7r + -7r~ pairs 
within a region of a few MeV (~ 50%<5m (7 ) around the peak. Not all of these pairs 
correspond to critical sigmas, but, according to the decay scenario described above, 
most of the critical sigmas certainly belong to this set. 

• Adding up the momenta of the selected charged pions form the momenta of the 
corresponding (fictitious) sigmas p a = p + + p- . 

• Perform a factorial moment analysis to the momenta of the fictitious sigmas. Due 
to the presence of the critical sigmas in this set, a strong intermittency effect shows 
up in the sigma sector (see Fig. 5, full circles) with the expected indices. 

7. Numerical results 

Having developed both a Monte-Carlo generator for critical events and an algorithm 
for the restoration of the critical sigma sector in an event-by-event analysis we shall now 
proceed to present our numerical results, showing the phenomenological impact of our 
theoretical investigations. Let us first specify values for our input parameters, which 
correspond to the critical exponent 5, the effective coupling g, the transverse radius R±, 
the proper time scale r, the size in rapidity A and the critical temperature T c . The first two 
parameters are fixed, on the basis of universality class arguments, to the values of the 3-D 
Ising model: 5 ~ 5 and g ~ 2 p2]JT3[] . The critical temperature is taken T c « 140 MeV. 
The size in rapidity is determined so as to meet the conditions at the RHIC collider: 
A « 11. Finally we choose: R± ~ 22 fm and r « 17 fm. It follows that the presented 
results are predictions for the event-characteristics in RHIC, provided that the critical 
point is there accessible. The above choice of R± and g leads us to the value R u ~ 13 fm. 
This value implies the formation of 2 clusters in the transverse momentum projection. 
Taking into account that A u « 0.55 we find 10 clusters in the rapidity projection. Thus 
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approximately 20 critical clusters are expected to be formed in RHIC conditions, adopting 
the cylindrical evolution scenario. Our numerical results are exhibited in a series of figures 
and pertain to a large set of critical events produced through the CMC algorithm and 
using parameter values appropriate for RHIC. We investigate separately the two possible 
geometries, cylindrical and spherical, for the evolution of the critical fluid. 

Case of cylindrical symmetry: We produce 10000 "critical" events and perform an 
event-by-event factorial moment analysis for the sigma clusters. Our results are depicted 
in Figs. 1-3. The mean multiplicity of sigmas turns out to be ~ 90. In Fig. la we present 
the inclusive rapidity spectrum of the sigmas for the set of the 10000 events. The uniform 
central plateau indicates the cancelling of the fluctuations in the inclusive distribution. 
In Fig. lb we show the inclusive distribution of the sigmas in px- 

In Figs. 2a,b results of the event-by-event factorial moment analysis in rapidity space 
and in transverse momentum space for the whole set of the 10000 "critical" events are 
displayed. In Fig. 2a we present the distribution of the slopes in a linear fit for ln(F 2 (1) (M)) 
as a function of ln(M) (M beeing the number of bins in rapidity) for all events. We 
observe a clear maximum near s 2 ~ 0.27 and a mean variance Ss^ ~ 0.05. This picture 
is consistent with the formation, in most events, of a complex structure with fractal 
dimension dp ~ 0.73 which is close to the mass dimension dp = | obtained for a single 
sigma cluster (see eq.([L5|)) in rapidity. In Fig. 2b the analogous analysis is made for the 
distribution of the sigmas in transverse momentum space. The corresponding histogram 

(2) 

gives the distribution of s 2 for the slopes determined after a linear fit to the logarithm of 
the two dimensional moment F 2 (M) as a function of the logarithm of M. We see, once 

(2) 

more, a peaked distribution with a maximum at about s 2 ~ 1.23. This distribution is 
slightly broader than the corresponding distribution in rapidity. The mean value (1.23) 

(2) 

indicates a fractal structure with dimension d F ~ 0.77 close to the mass dimension of a 
single cluster in p? (dp = |). For completeness we mention that for the azimuth angle in 
transverse momentum space we get a totally uniform inclusive distribution, as one would 
expect. 

Putting pions into the picture we give, in Figs. 3a,b, the results of the event-by-event 
factorial moment analysis (a) in rapidity and (b) in px of the negative pion momentum 
distribution. To determine the decay of sigmas into pions according to the distribution 
P(Mo) we use the parameter values M a « 285 MeV and 8m„ ~ 5 MeV. As pion mass 
we have used = 140 MeV. For the mean value of s 2 i n the rapidity space analysis of 
7i we find < 4 1} >^ 0.07 (see Fig. 3a) while for the case of the transverse momentum 
we find < 4 2) >~ 0.17 (see Fig. 3b). These values should be compared with the values 
< s% >~ 0.27 and < >~ 1.23 respectively, found in the event-by-event analysis 
of the sigmas. We see a clear tendency for suppression of the dynamical fluctuations in 
favour of the statistical ones in the pionic sector. 

To recover the dynamical fluctuations one has to apply the sigma-recontruction al- 
gorithm described in the previous paragraph. To investigate the suppression and the 
restoration of the critical fluctuations we compute the dependence of < s 2 > f° r the 
reconstructed sigmas on the quantity 8m a , keeping a fixed value M a = 290 MeV. The 
analysis is performed for classes (constant 5m a ) with 500 events each. Our results are pre- 
sented in Fig. 4 (full squares). In the same plot we show the best fit using the Boltzmann 
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function f(x) = 1+ ^ ( 1 3 ._f o 2 )/d + A 2 (dashed line). The fit parameters converge to the values: 
Ai = 2.12, A 2 = 0.35, x = 0.08 and d = 1.1. The nonvanishing value of A 2 shows that, 
even if 5m a becomes very large, a clear intermittency signature (< s 2 > > 0.35) is 
present in the reconstructed events. Of course the asymptotic value of < s 2 ^ > depends 
on M a and decreases rapidly as M a is increased. 

To illustrate this effect more transparently we present in Fig. 5 the second factorial 
moment in transverse momentum space for a typical critical event taking M a = 290 MeV, 
bm a = 5 MeV and m w = 140 MeV. We see that although the fluctuations in n~ (open 
up triangles) are strongly suppressed they are practically revealed in the sector of the 
reconstructed sigmas (full circles). 

Spherically symmetric evolution: In Figs. 6a,b we give the rapidity (a) and transverse 
momentum (b) inclusive distributions for the sigmas in a set of 10000 Monte-Carlo gen- 
erated events. The distribution in rapidity is very interesting. The formation of a single 
source (cluster) of sigmas in the collision process leads to an almost Gaussian shape which 
is very clearly distinguished from the plateau emerging in the case of many sources (cylin- 
drical symmetry). Such differences in the rapidity distribution are also observed between 
p — p and heavy-ion collision experiments far from the critical point. Our explanation is 
that the multiparticle production in p — p is dominated by the formation of a single source 
while in heavy-ion collisions there are many sources involved (supporting the scenario of 
longitudinal expansion). 

Finally, in Figs. 7a,b we present the factorial moment analysis in the sigma sector for 

the set of the 10000 MC events presented in Figs. 6a,b. We find < s 2 >= 0.14 for 

(2) 

the rapidity distribution and < s 2 >= 0.99 for the transverse momentum distribution. 
A suppression of the fluctuations in the rapidity space is clearly observed. This is a 
peciularity of the spherical geometry as discussed in the previous section. 

In summary we have shown that even if the sigmas, which are expected to be formed 
in a heavy ion collision experiment, as the critical point is approached, decay well above 
the two-pion threshold, it is possible, by investigating the momentum distribution of the 
produced charged pions, to (partially) reconstruct the critical sigma sector and to restore 
the corresponding critical fluctuations. 

8. Concluding remarks 

In this work we have explored the idea that pion production in heavy ion collisions 
may have a component whose source can be asrcibed to a critical QCD phase. Our 
working hypothesis has been that density fluctuations near the QCD critical point should 
leave their imprints on the observed pion spectrum. Encouraged from the accumulating 
evidence, through microscopically-based theoretical investigations, lending support to a 
phase diagram in which critical QCD behaviour does occur we focused our efforts on 
establishing a linkage between critical fluctuations and intermittency patterns in pion 
production. Such a proposal pertains, of course, to one of two components characterizing 
the pion spectrum in A + A collisions, as conventional production mechanisms are also 
expected to take place. 

The methodology employed in this work centered itself on the statistical mechanical as- 
pects of the problem. We have chosen to rely on the strategy of critical cluster formation 
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which seems to be the most advantageous one given the finite size aspects of the system 
under study. In the context of the currently accepted theoretical proposal that the uni- 
versality class of the critical QCD system is represented by the 3-D Ising model, we have 
both established connections with pion production ascribed to the aforementioned un- 
convential component and explored its observability in terms of specific filters associated 
with a reasonably wide window for the a-mass. Employing Monte-Carlo type numeri- 
cal methodologies we have arrived at specific testable predictions which can be checked 
against forthcoming data in heavy ion experiments. In practical terms our present under- 
taking can be viewed as having a dual prospective. On the one hand, it aims at isolating a 
specific type of signal in the constitution of the final states with an unconventional origin 
which would give evidence in favor of an antecedent QGP phase. On the other, it could 
serve as a guide towards the identification of the source (critical point?) responsible for 
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The fact that the pions constitute the focal point of our approach is also significant 
from a purely phenomenological point of view addressing itself directly to the study of the 
constitution of the multiparticle spectrum (particle multiplicities and particle ratios) pro- 
duced in very high energy collisions, especially A In this connection, thermodynamical 
analyses from the hadronic perspective which base themselves on thermal equilibration 
assumptions for the produced multiparticle system are of special interest. Pioneering 



studies in this direction [20] have assumed an ideal hadron gas model and have produced 
satisfactory results, lending credit to the thermal equilibration premise, but are in no 
position to evaluate whether or not the source of the resulting multiparticle system is 
located inside or outside (QGP?) the hadronic domain. 

Recently, the same type of analysis has been conducted [2~I| within the framework of 
Hagedorn's Statistical Bootstrap Model (SBM) |22| once the latter was extended so as 
to include the strangeness quantum number. This approach not only incorporates, via 
a self-consistent logic, interactions into the resulting thermodynamical system but also 
sets limit for the hadronic phase of matter. According to the results of there is 
an accumulation, in the plane of thermodynamical variables specified by temperature 
and light quark chemical potential, of points from S+S, S+Ag and (preliminary) Pb+Pb 
experimental data from SPS at CERN on either side of the curve that separates the 
hadronic world from a different phase of matter. Matching such results from the QCD 
side certainly specifies a task worthwhile pursuing (suggesting a first order transition). 

Looking ahead, we assess as the most pressing question in relation to our proposed 
scheme to be the following. If, according to the present theoretical evidence, the second 
order phase transition achieved by the QGP system in the temperature-baryon chemical 
potential plane is a critical point marking the end of a first order transition curve, are there 
realistic prospects of observing the component contributing to pion production presently 
investigated? Clearly, the answer depends as to whether one can define a reasonably 'wide' 
window around such a critical point within which visible implications can be formulated. 
This constitutes a problem of interest which is currently being investigated. At the same 
time, we believe that the methodology developed in this paper has a wider impact as 
it can be easily adopted to any future scenario concerning the, still open, issue of QCD 
universality classes implicated as being realized in A + A collisions. 
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Figure 1. (a) The rapidity distribution of the a-particles in 10000 Monte-Carlo (MC) 
generated critical events (cylindrical geometry is assumed to be valid), (b) The corre- 
sponding inclusive transverse momentum distribution. 
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Figure 2. (a) The distribution of the slope s 2 resulting from a linear fit of ln(F 2 ) 
vs. ln(M) in an event-by-event factorial moment analysis of the a-particles in rapidity 
for the 10000 MC events (M is the number of bins), (b) The distribution of the slope 

(2) 

s 2 obtained in the analogous analysis in transverse momentum space for the same set of 
events as in Fig. 2a. 
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Figure 3. (a) The distribution of the slope resulting from a linear fit of li^i*^) vs. 
ln(M) in an event-by-event factorial moment analysis of the negative pions in rapidity 

(2) 

for the 10000 MC events (the same as in Figs. 1,2). (b) The distribution of the slope s 2 

(2) 

resulting from a linear fit of ln(F 2 ) vs. ln(M) in a 2-D event-by-event factorial moment 
analysis of the negative pions in transverse momentum space. 
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Figure 4. The dependence of < s 2 > on Jm^- (full squares) calculated with sets each 
consisting of 500 critical events generated with the cylindrically symmetric algorithm. Mo- 
is taken 290 MeV. A fit with the Boltzmann function (dashed line) is also shown. 
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Figure 5. The second factorial moment in transverse momentum space for a critical 
event (negative pions) obtained through the cylindrically symmetric CMC generator and 
corresponding to M a = 290 MeV, 5m a = 5 MeV (open up triangles). The corresponding 
moment in the sigma-sector obtained through the above described reconstruction algo- 
rithm is also shown (full circles). The restoration of the critical fluctuations is clearly 
displayed. 
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Figure 6. Inclusive distributions in (a) rapidity and (b) magnitude of the transverse 
momentum of the a-particles for 10000 Monte-Carlo generated critical events assuming 
spherical symmetry. 
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Figure 7. (a) The distribution of the slope sj> resulting from a linear fit of ln(F 2 ) vs. 
ln(M) in an event-by-event factorial moment analysis of the sigmas in rapidity for 10000 
MC events with spherical evolution, (b) The distribution of the slope resulting from 
a linear fit of lnf^ 2 "*) vs. ln(M) in a 2-D event-by-event factorial moment analysis of the 
sigmas in transverse momentum space for the same set of events as in (a). 



